Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 November 2012 (MN M£X style file v2.2) 



Simulating the assembly of galaxies at redshifts z = 6 — 12 



Pratika Dayal 1 *, James S Dunlop 1 , Umberto Maio 2 & Benedetta Ciardi 3 

' ^ 1 SUPA], Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh, EH9 3HJ, UK 
' _ ^ l 2 Max- Planck- Institut fur extraterrestrische Physik Giessenbachstrasse 1, 85748 Garching, Germany 
■ 3 Max Planck Institut fur Astrophysik, Kari-Schwarzschild-Strasse 1, 85741 Garching, Germany 



ABSTRACT 

We use state-of-the-art simulations to explore the physical evolution of galaxies in 
the first billion years of cosmic time. First, we demonstrate that our model, without any 
tuning, reproduces the basic statistical properties of the observed Lyman-break galaxy 
(LBG) population at z = 6 — 8, including the evolving ultra-violet (UV) luminosity 
*-£h ■ function (LF), the stellar- mass density (SMD), and the average specific star- formation 

rates (sSFR) of LBGs with Muv < —18 (AB mag). Encouraged by this success we 
present predictions for the behaviour of fainter LBGs extending down to Muv — —15 
(as will be probed with the James Webb Space Telescope) and have interrogated our 
simulations to try to gain insight into the physical drivers of the observed population 
evolution. We find that mass growth due to star formation in the mass-dominant 
progenitor builds up about 90% of the total z ~ 6 LBG stellar mass, dominating 
over the mass contributed by merging throughout this era. Our simulation suggests 
that the apparent "luminosity evolution" depends on the luminosity range probed: the 
steady brightening of the bright end of the LF is driven primarily by genuine physical 
luminosity evolution and arises due to a fairly steady increase in the UV luminosity 
(and hence star-formation rates) in the most massive LBGs; e.g. the progenitors of the 
z ~ 6 galaxies with Muv < —18.5 comprised ~ 90% of the galaxies with Muv < —18 
at z ~ 7, and ~ 75% at z ~ 8. However, at fainter luminosities the situation is more 
complex, due in part to the more stochastic star-formation histories of lower-mass 
objects; the progenitors of a significant fraction of z ~ 6 LBGs with Muv > — 18 
were in fact brighter at z ~ 7 (and even at z ~ 8) despite obviously being less 
massive at earlier times. At this end, the evolution of the UV LF involves a mix of 
positive and negative luminosity evolution (as low-mass galaxies temporarily brighten 
then fade) coupled with both positive and negative density evolution (as new low- 
mass galaxies form, and other low- mass galaxies are consumed by merging). We also 
predict the average sSFR of LBGs should rise from sSFR ~ 4.5 Gyr -1 at z ~ 6 to 
sSFR - 11 Gyr" 1 by z ~ 9. 

Key words: Galaxies: evolution - high-redshift - luminosity function, mass function 
- stellar content 



1 INTRODUCTION 

In the standard Lambda Cold Dark Matter (ACDM) sce- 
nario, the first structures to form in the Universe were low- 
mass dark-matter (DM) halos, and these building blocks 
merged hi erarchically to form su ccessively larger struc- 
tures (e.g. iBlumenthal et all 1 1984 ). While the initial con- 
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ditions (at z fa 1100) for the formation of such structures 
are well established as being nearly scale-invariant, small- 
amplitude density fluctuations (e.g. iKomatsu et al.ll2009T ). 
until recently there has been relatively little direct obser- 
vational information on the emergence of the first galax- 
ies. These early galaxies changed the state of the inter- 
galactic medium (IGM) from which they formed: they pol- 
luted it with heavy elements and heated and (re)ionized it, 
thereby affecting the evolution of all subsequent generations 
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of galaxies dBarkana fc Loebll20oil ; ICiardi fc Ferraral 120051 ; 
iRobertson et al.ll201Ch ~ 

As recently reviewed by iDunlgpl |2012f) . until the ad- 
vent of sensitive near-infrared imaging with the installation 
of Wide Field Camera 3 (WFC3/IR) on the Hubble Space 
Telescope (HST) in 2009, there existed at most one con- 
vinci ng galaxy candidate at redshifts z > 7 (Bo uwens et alj 
12004 ) . Now, however, over one hundred galaxies have been 
uncovered with HST in the redshift range 6.5 < z < 8.5 (e.g . 
Qesch et al.ll2010|;lBouwens et al.ll2010ai;lMcLure et al.l2010l; 
Finke lstein et al.ll2010l: iMcLure et al.112011 Lorenzoni et al.l 
201ll ; lBouwens et al.ll201ll ; lOesch et al.ll2012h . and the new 



generation of ground- based wide-field near-inf rared surveys, 
such as Ultra VISTA (|McCracken et alj|2012h are reaching 
th e depths required to reveal more luminous galaxies at z ~ 
7 dCastellano et al.|[201Qbl ; lOuchi et al. I l20ld iBowler" et al.l 
12012ft . 

The assembly of significant samples of galaxies within 
< 1 Gyr of the big bang has enabled the first mean- 
ingful determinations of the form and evolution of the 
UV LF of galaxies over the redshift range 6 < z < 8, 
reaching down to absolute magnitudes Muv — —18 



(IMcLure et al.l l20ld; 


Oesch et al. 


l20ld; iFinkelstein et al.1 


20ld; Bouwens et al. 


2011; Finkelstein et al. 2012a; 



Bradley et all l2012f ). It has also facilitated the study 
of the typical properties of these young galaxies , both 
through follow-u p spec troscopy l|Schenker et al.l |2012| ; 
ICurtis-Lake et ahl l2012bf ), and from analysis of their 
average broad-band (H ST + Spitzer) spec t ral energy dis- 



2010b; Gonzalez et al.1 


2010; 


Finkelstein et al. 




2010]; 


Wilkins et al.l l201ll; IMcLure et al.l l201ll; iDunloo 


et al.l 


20121; Bouwens et al.1 2012; Finkelstein et al.1 


2012b: 


Curtis-Lake et al.1 2012a; 


IStark et al. |2012|; Labbe et al.1 


2012|). 



The results of these studies are still a matter of vigorous 
debate, and will undoubtedly be refined substantially over 
the next few years, as the HST/ Spitzer and ground-based 
near-infrared imaging database expands and improves, and 
deep multi-object near-infrared spectroscopy becomes rou- 
tine. Nevertheless, the current situation can be broadly sum- 
marized as follows. First, there is good agreement over the 
basic form and evolution of the galaxy UV LF at z ~ 6, 7, 8, 
albeit the current data do not allow the degeneracy between 
the Schechter function parameters to be unambiguously bro- 
ken. Second, there is as yet no meaningful information on the 
number density of galaxies at z > 8.5. Third, there is some 
tentative evidence for a drop in the prevalence of Lyman-a 
line emission at z ~ 7, as arguably expected from an in- 
creasingly neutral ISM. Fourth, various studies suggest that 
the specific star formation rate (sSFR — star-formation 
rate/stellar mass) of star- forming galaxies, after rising from 
z = to z ~ 3, then remains remarkably constant out to 
z ~ 7. Fifth, there is growing evidence that this last result 
may need to be modified in the light of increasingly strong 
contributions from nebular line emission with increasing red- 
shift. Sixth, the UV continua of the highest redshift galaxies 
is fairly blue (indicative of less dust than at lower redshifts) , 



but there is as yet no evidence for exotic stellar populations, 
as might be revealed by extreme UV slopes. 

Despite the remaining uncertainties, it is clear that 
these new observational results already provide an excel- 
lent test-bed for the increasingly complex galaxy formation 
and evolution models/simulations tha t have recently been 
developed in a number of studies (e.s;. [Finlator et al.ll2007l; 



Daval et al.l 120091. 120101; I Salvatcrr a et al 



2011 



Daval et al.l 



20111; iForero-Romero et alj l201ll; IDaval fc Ferraral I2OI2T: 



Yaiima et a"i1l2012l '). 



In this study, our aim is not only to test the ability 
of the latest galaxy formation models to reproduce these 
new observational results at high-redshift, but also to use 
our simulations to obtain a better physical understanding of 
how the emerging galaxies evolve in luminosity and grow in 
stellar mass. We choose z ~ 6 as the baseline, and explore 
the evolution and physical properties of the progenitors of 
z ~ 6 galaxies out to z ~ 12. We explore the properties 
of the currently observable galaxies with Muv > —18, but 
also present predictions reaching over an order-of-magnitude 
fainter, to Muv — —15. Such faint galaxies are of interest be- 
cause i) they are extremely numerous, and thus provide good 
number statistics for observational predictions, ii) they are 
expected to have low dust content (IDaval fc Ferraral |2012| . 
see), simplifying the predictions for rest-frame UV observa- 
tions, iii) current evidence suggest it is these faint galax- 
ies which likely provide the b ulk of the photons responsi- 
ble for Hydrogen reio nization l|Choudhurv fc Ferraral [20071 : 
ISalvaterra et al.ll201l[ ). iv) such galaxies are within reach of 
planned observations with the James Webb Space Telescope 
(JWST) at the end of the decade, and v) the simulation of 
such faint dwarf galaxies requires exquisite mass-resolution, 
achievable with the simulation utilised in this work. Specif- 
ically, we use state-of-the-art simulations with a box size of 
lOcMpc (comoving Mpc), and include molecular hydrogen 
cooling, a careful treatment of metal enrichment and of the 
transition from Pop III to Pop II stars, along with modelling 
of supernova feedback, as explained in Sec. [2] that follows. 



2 COSMOLOGICAL SIMULATIONS 

The simulation used in this work is no w briefly des c ribed , 
and i ntere sted readers are referred to iMaio et al.l (|2007l . 
20091, l201Ch and ICampisi et all l|201ll ) for complete details. 
We use a smoothed particle hydrodynamic (SPH) simula- 
tion carried ou t using the TreePM-SPH code GADGET-2 
(Springcl 2005). The cosmological model corresponds to the 
ACDM Universe with DM, dark energy and baryonic density 
parameter values of (f2 m , Qa, fib) = (0.3, 0.7, 0.04), a Hubble 
constant Ho — WOh — 70kms _1 Mpc _1 , a primordial spec- 
tral index n s — 1 and a spectral normalisation erg = 0.9. The 
periodic simulation box has a comoving size of 10/i _1 cMpc. 
It contains 320 s DM particles and, initially, an equal number 
of gas particles. The masses of the gas and DM particles are 

3 x 10 5 /t _1 M© and 2 x 10 6 /i _1 M Q , respectively. The max- 
imum physical softening length is set to 0.5/i -1 Kpc and is 
allowed to drop (at most) to half of the gravitational soft- 
ening. 
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The code includes the cosmological evolution of e 



H, ET, He, He + , He" 1 



H 2 , Hj-, D, D + , HD, HeH 



ijYoshida et al.l 120031 : iMaio et al.l l2007l . 120091 1 . PopIII and 
PopII/I star formation ac cording; to the corresp onding ini- 



tial mass function (IMF; iTornatore et al 



2007), _gas cool- 



ing from resonant an d fine-structure lines ( Maio et alj|2007f ) 
and feedback effects (|Springel fc Hernquistll2003l ). The runs 
track individual heavy elements (e.g. C, O, Si, Fe, Mg, S), 
and the transition from the metal-free PopIII to the metal- 
enriched PopII/I regime is determined by the underlying 
metallicity of the mediu m, Z, compared with the critical 
" 4 Z(7, (iBromm fc Loebll2003l, and refer- 
Schneider et ali (|2003, 2006|) who 



value of Z cr it = 10" n & 
ences therein) , but see also \ 
find a much lower value of Z cr i t = 
precise value of Z cr i t is not crucial; 



10"° Z . However, the 
the large metal yields 
of the first stars easily pollute the medium to metalli city 
values abov e ~ 10~ 4 - 10~ 3 Z l|Maio et al.ll2010l . l201ll ). If 
Z < Z cr it, a lSalpeter] (|l955h IMF is used, with a mass range 
100 — 500 M©; otherwise, a standard Salpeter IMF is used 
in the mass ra nge 0.1 — 100 Mq, and a SNII range is used 
for 8 - 40 M (IBromm et al.ll2009l ; IMaio et al.ll201ll ). 

The chemical model follows the detailed stellar evolu- 
tion of each SPH particle. At every timestep, the abun- 
dances of different species are consistently deriv e d us - 
ing the lifetime function (Pad ovani fc Matteuccil 1 19931 ) 
and metallicity-dependant stellar yields; the yields from 
SNII, AGB stars, SNIa and pair instability SN (PISN; 
for primordial SN n ucleos y nthesis) have been taken from 
IWooslev fc Weaver! d 19951) . van den Hoek fc Gro cnewegcn 
'l997 h , iThielemann et"ai1 l|2003l ) and iHeger fc Wooslev 



2002), respectively. Metal mixing is mimicked by smooth- 
ing metallicities over the SPH kernel. Pollution is driven 
by wind feedback, w hich causes metal spreading over ~ kpc 
scales at each epoch (|Maio et al.ll201lh . 

Galaxies are recognized as gravitationally-bound groups 
of at least 32 total (DM+gas+star) particles by running a 
friends-of-friends (FOF) algorithm, with a comoving linking- 
length of 0.2 in units of the mean particle separation. Sub- 
structures are identified by using th e SubFind algorithm 
l|Springel et al.ll200ll ; iDolag et al.ll2009T l which discriminates 
between bound and non-bound particles. Of the galaxies 
identified in the simulation, at any redshift, we only use 
'well-resolved' galaxies in our calculations, i.e. those that 
have at least 10 star particles; this corresponds to galaxies 
with more than 145 total particles at any of the redshifts 
studied, 2 ~ 6 — 12. For each such well-resolved galaxy, we 
obtain the properties of all its star particles, including the 
redshift of, and mass/metallicity at formation; we also com- 
pute global properties including the total stellar/gas/DM 
mass, mass-weighted stellar/gas metallicity and the instan- 
taneous star-formation rate (SFR). 



2.1 Identifying LBGs 

To identify the simulated galaxies that could be observed 
as LBGs, we start by computing the UV luminosity (or ab- 
solute magnitude) of each galaxy in the simulation snap- 
shots at z ~ 6 — 12, roughly spaced as Az ~ 1. We con- 
sider each star particle to form in a burst, after which it 



evolves passively. As explained above, the simulation used 
in this work has a detailed treatment of chemodynamics, 
and consequently we are able to follow the transition from 
metal- free PopIII to metal-enriched PopII/I star formation. 
If, at the time of formation of a star particle, the metallic- 
ity of its parent gas particle is less than the critical value 
of Z cr it = 10 -4 Zq, we use the PopIII SED, taking into ac- 
count the mass of the star particle so formed, and assuming 
that the emission properties of PopIII stars remain constant 
during their re latively short lifetime of about 2.5 x 10 6 yr 
|Schaererll2002T ). On the other hand, if a star particle forms 
out of a metal-enriched gas particle (Z > Z cr it), the SED is 
computed via the po pulation synthesis code STARBURST99 
l|Leitherer et al.lll999l) , using its mass (M»), stellar metallic- 
ity (Z*) and age (t*); the age is computed as the time differ- 
ence between the redshift of the snapshot and the redshift 
of formation of the star particle. Although the luminosity 
from PopIII stars has been included in our calculations, it 
is negligibl e compared to that co ntributed by PopII/I stars 
(see Fig. 9. ISalvaterra et al1l201ll ). The composite spectrum 
for each galaxy is then calculated by summing the SEDs 
of all its (metal-free and metal-enriched) star particles. The 
intrinsic continuum luminosity, L c n ', is calculated at rest- 
frame wavelengths A = 1350, 1500, 1700 A at z ~ 6, 7, 8 
respectively; though these wavelengths have been chosen for 
consistency with J-band observations, using slightly differ- 
ent values would not affect the results in any significant way, 
given the flatness of the intrinsic spectrum in this relatively 
short wavelength range. 

Dust is expected to have only a minor impact on the 
faint end of the UV LF at z > 6 where {E(B - V)) < 0.05, 
as a result of the low stellar mass, age and metall icity of 
these early galaxies (see also ISalvaterra et al.l 1201 lh . Nev- 
ertheless, we still self-consistently compute the dust mass 
and attenuation for each galaxy in the simulation box; we 
note that while the metal enrichment and mixing have been 
calculated within the simulation as explained in Sec. O 
the dust masses and attenuation are calculated by post- 
processing the simulation outputs. Dust is produced both 
by SN and evolved stars in a galaxy. However, the contri- 
bution of AGB stars becomes progressively less important 
towards very high redshifts, because the typical evolutionary 
time-scale of these stars (Jj 1 Gyr) becomes longer than the 
age of the Universe above 2: ^ 5 .7 ( Todini fc Ferrarall200ll ; 
iDwek et al.1 120071 : IValiante et al.ll2009l ). We therefore make 
the hypothesis that the dust present in galaxies at z ^ 6 is 
produced solely by SNII. The total dust mass present in each 
galaxy is then c omputed assuming: (i) 0.4 M^ of dust is pro- 
duced per SNII llTodini fc Ferrarail200ll : iNozawa et aHl2007l : 
iBianchi fc Schneiderl 120071 ). (ii) SNII destro y dust in for- 
ward shocks with an efficiency of about 20% l|McKeel ll989: 
ISeab fc Shuli| [T983). (iii) a homogeneous mixture of gas and 
dust is assimilated into further star formation (astration), 
and (iv) a homogeneous mixture of gas and dust is lost in 
SNII powered outflows. 

To transform the total dust mass into an optical depth 
to UV continuum photons, we assume the dust to be made 
up of carbonaceous grains and spati ally distributed in a 
screen of radius = r g — 4.5Ar2oo (jFerrara et al.1 120001 ) ■ 
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Figure 1. The UV LFs of galaxies at z ~ 6,7, 8 as marked in each panel. In all panels, histograms show our theoretical results for 
galaxies with Mjjy < —15, with error bars sho wing the associated p oissonian errors, a nd points show the o bserved data. The observed 
LBG UV LFs have been taken fro m: (a) z ~ 6: | Bouwens et al] ||2007| , fille d circles') and | McLure et aT] ||2009| . fi lled triangles); (b ) z ~ 7 : 
lOesch et al.l ||2010 | . filled squares), IBouwens et alj ll2010bl , empty circles ) . IBouwens et alj J201 ll . filled cir cles). |Castellano et alj J2010tj . 
empty t riangles) . IvicLure et al"] ||2010| , filled triangles) and |Bow ler et al] ||2012| , empty squares); (c) z ~ 8: |Bouweiis**eT*aLTl20*10bl . empty 
circles). IBouwens et alj ll201ll . filled circles). iMcLure et al] I2OIOI . filled triangles) and iBradlev et aLI ll2012l , empty squares). The vertical 
short (long) dashed lines in each panel show approximate effective detection limits for HST (JWST) near-infrared imaging. 



where r a is the radius of the gas distribution, A = 0.05 is 
the spin parameter, and r2oo is the virial radius, assuming 
that the collapsed region has an over-density which is 200 
times the critical d ensity at the reds hift considered. The 
reader is referred to lDaval et al.l l|2010l ) for complete details 
of this calculation. The observed UV luminosity can then 
be expressed as L° bs = L* ni x / c , where f c is the fraction of 
continuum photons that escape the galaxy, unattenuated by 
dust; UV photons (A > 1500A) are unaffected by the IGM, 
and all the UV photons that escape a galaxy unattenuated 
by dust reach the observer. 



3 PREDICTED OBSERVABLES FOR LBGS 

Once the UV luminosity of each galaxy in the simula- 
tion snapshots at z ~ 6 — 12 has been calculated as ex- 
plained above, we can study the galaxies that would be de- 
tectable as LBGs. Although the current near-infrared ob- 
servable limit of the deepest HST surveys corresponds to 
Muv (-17.75, -18.0, -18.25) at z ~ (6, 7, 8), as a result 
of the exquisite simulation resolution we are able to study 
LBGs that are an order of magnitude fainter (Muv < —15); 
these faint galaxies will be detectable with future facilities 
such as the JWST. As a first test of the simulation, we now 
compare the UV LFs, specific star-formation rate (sSFR) 
and the stellar mass density (SMD) of the simulated LBGs 
to the observed values. 



3.1 UV Luminosity Functions 

We start by building the UV LFs for all galaxies with 
Muv < —15 in the simulation boxes at z ~ 6, 7 and 8. 
It is clear that the luminosity range sampled by the obser- 
vations and the theoretical model overlap only in the range 
Muv ^ -17.75 to -19.5, as shown in Fig. □ This is be- 
cause the data are limited at the faint end by the deep- 



est HST+WFC3/IR data obtained prior to the upcoming 
UDF12 campaign (HST GO 12498, PI: Ellis), while the 
models are limited at the bright end by the size of the simu- 
lated volume needed to achieve the exquisite mass resolution 
(~ 1O 5 M0 in baryons) required to model the PopIII to II 
transition in high-z galaxies. Nevertheless, in the magnitude 
range of overlap, both the amplitude and the slope of the 
observed and simulated LFs are in excellent agreement at 
all three redshifts. This is a notable success of the model, 
given that once the SEDs and UV dust attenuation of each 
galaxy in the simulated box have been calculated as de- 
scribed above, we have no free-parameters left to vary when 
computing the predicted UV LF. 

The simulated UV LFs reproduce well the observed pro- 
gressive shift of the galaxy population towards fainter lu- 
minosities/magnitude (M^v) and/or number densities (</>*) 
with increasing redshift from z ~ 6 to 8. Observationally 
the physical driver of this evolution is not yet clear, with the 
Schechter fun ction fits of some authors favouring luminos- 
ity evolution ( Mannucci ct al. 2007; Castcll ano et al]|2010al : 
iBunker et al] 2010l ; Bouwens et al]|201ll ). while the results 



of other studies appear to favour prima rily density evolution 
|Ouchi et al]|2009l ; lMcLure et alj201Ch . However, by tracing 
the histories of the individual galaxies in the our simulation 
we can hope to shed light on the main physical drivers of 
this population evolution (see Sec. 5). 

Our simulations also reproduce well the observed steep 
faint-end slope of the UV LF, consistent with a ~ —2 at 
all three redshifts. There is, however, still considerable de- 



(lOesch et al]|2010l; 


McLure et al]|2010l; IBouwens et al]|201ll; 


IBradlev et al.ll2012 


). The faint-end slope of the LF has im- 



portant implications for reionization since these faint galax- 
ies are expect ed to be the major sources of H 1 ionizing pho- 
tons (see e.g. IChoudhurv fc Ferraral 120071; iRobertson et alj 
l201fj|Salvaterra" et al ]|201ll ; lFerrara fc Loebll2012l ). The cur- 
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Figure 2. Stellar mass density (SMD) as a function of red- 
shift. Filled squares (triangles) show our theoretical predictions 
for galaxies with Muv < —15 (Muv < —18) at each redshift, 
with error bars showing the associated poissonian errors. Empty 
points show the magnitude-limited SMD values inferred observ a- 
tionally for galaxies with M uv < -18 by lLabbe et all fcOlOal lbl. 
empty stars) and IGonzalez et alj d201l|. empty tr i angles ) ; empty 
squares show the SMD inferred by IStark et all J2012h for the 
same magnitude limit after correcting the stellar masses for neb- 
ular emission-line contributions to the broad-band fluxes (assum- 
ing that the nebular line rest-frame equivalent width evolves with 
redshift). 



rent uncertainty in the value of a should be clarified some- 
what by UDF12 (Mclure et al., in preparation), and should 
be resolved by JWST with its forecast ability to reach ab- 
solute magnitude limits of Muv — (—15.25, —15.5, —15.75) 
at z ~ (6, 7, 8) respectively. 



3.2 Stellar mass density 

The stellar mass, M», is one of the most fundamental prop- 
erties of a high-z galaxy, encapsulating information about 
its entire star-formation history. However, achieving accu- 
rate estimates of M* from broad-band data is difficult be- 
cause it depends on the assumptions made regarding the 
IMF, SF history, strength of nebular emission, age, stellar 
metallicity and dust; the latter three parameters are degen- 
erate, adding to the complexity of the problem. Although 
properly constraining M* ideally requires rest-frame near 
infra-red data which will be provided by future instruments 
such as MIRI on the JWST, broad-band HST+Spitzer data 
have already been used to infer the contribution of galax- 
ies brighter than Muv = — 18 to the growth in total stellar 
mass density (SMD = stellar mass p er unit volume) wi th 

; lLabbe et ai1l2010al lbl; 



decreasing redshift 


Stark et alJl200Sl 


IGonzalez et al.ll2011 


). Encouragingly, 



observed growth in SMD is reproduced very well by inte- 



Figure 3. Specific star-formation rate (sSFR) as a function of 
redshift. Filled squares (triangles) show the theoretical results for 
galaxies with Muv < — 15 (Muv < —18) at each redshift, with 
error bars showing the associated poissonian errors. E mpty points 
show the sSFR value s inferred obser v ationa lly by IStark et alj 
||2009| , empty circles), IGonzalez et alj J 20 1 dl . empt y triangles); 
empty stars and squares show the sSFR inferred by Stark et al. 
|2013), after correcting the stellar masses for nebular emission 
lines assuming that the nebular line rest-frame equivalent width 
at z ~ 4 — 7 is the same as that as z ~ 3.8 — 5, and alterna- 
tively that the nebular line rest-frame equivalent width evolves 
with redshift, respectively. The solid and dashed lines show the 
evolution of sSFR as inferred from the simulation for galaxies 
with M, = 10 6_8 Mq and M* > 10 8 M Q , respectively. 



grating the stellar masses of the simulated galaxies brighter 
than Muv < —18; the theoretically calculated SMD of these 
galaxies drops to zero at 2 > 9 since there are no galaxies 
massive enough to be visible with this magnitude cut in the 
volume simulated. Further, as a result of their much larger 
numbers, galaxies with — 18 ^ Muv ^ — 15 contain about 
1.5 times the mass as compared to the larger and more lumi- 
nous galaxies that have been observed as of date; the JWST 
will be instrumental in shedding light on the properties of 
such faint galaxies, in which most of the stellar mass is locked 
up at these high-z. 



3.3 Specific star formation rates 

The specific star-formation rate (sSFR) is an important 
physical quantity that compares the current level of SF to 
the previous SF history of a galaxy. It also has the advan- 
tage of being relatively unaffected by the assumed IMF. 
There is now a considerable body of observational evidence 
indicating that the typical sSFR of star-forming galaxies 
rises by a factor of ab out 40 between 2 = and 2 ~ 2 
(e.g. iDaddi et al.ll2007l ). but then settles to a value con- 
sistent with 2 - 3Gyr _1 at 2 ~ 3 - 8 (IStark et all 120091 ; 
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Table 1. For the rcdshift, z, given in Column 1, we show the following ratios: SM^mbW/SM,^^ (column 2), 
£M» iMB (z)/£M» iaiipro9 (z) (column 3), EM, |O j Ipr0 g(z)/SAf» lS!= 6 (column 4), SM sf (z,z - l)/SAf,, e=6 (column 5), T,M acc (z,z - 
l)/SAf* iZ — e (column 6) and EM KC (z,2 — l)/T,M s f(z, z — 1) (column 7). Here, SM* >Z — 6 represents the total stellar mass at z ~ 6, 
summed over all z ~ 6 LBGs in the simulation (and has a value EM, jZ= s = 5.02 X 10 10 Mq), SM, mb( z ) is the stellar mass summed 
over all the major branch progenitors of z ~ 6 LBGs at the given redshift, EM* a jj pro „(«) is the total stellar mass in all the other 
progenitors of z ~ 6 LBGs (i.e. excluding the major branch one) at the redshift z, and SM s t(z, z — 1) and T,M acc (z, z — 1) represent 
the total stellar mass assembled by the major branch progenitors of z ~ 6 LBGs by star formation and mergers respectively, between 
the redshifts z and z — 1 (see Sec. 14.31 for details). 
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Gonzalez et~al]|2010l : iMcLure et afl|201ll ; IStark et al.ll2012i ; 
Labbe et al.ll2012l ). 



This constancy of sSFR at high redshift has proved 
somewhat unexpected and difficult to understand, given 
that theoretical models predict that sSFR should trace 
the baryonic infall rate which scales as (1 + z) 2,25 
l|Neistein fc Dekell lioosl ; IWeinmann et all 1201 ll ); according 
to this calculation, the typical sSFR should increase by a 
factor of abou t 9 over the redshi ft ran ge z ~ 2 — 7 . How - 
ever, recently, IStark et al.1 (|2012h and iLabbe et all (|2012T l 
have suggested that at least part of this discrepancy might 
be due to a past failure to account fully for the (poten- 
tially growing) level of nebular emission-line contamina- 
tion of broad-band photom etry at h igh-z, as su ggested by 



ISchaerer fc de Barrosl (|2009l ); indeed. IStark et aU (|2012ft es- 
timate that the emission-line corrected sSFR may rise by a 
factor of about 5 between z ~ 2 — 3 and z ~ 7. 

From the simulation snapshots, we find that typical 
sSFR decreases with increasing M*. However, this mass 
dependence is relatively modest, with the sSFR of galax- 
ies with Mi. — 10 6_8 Mq being only ~ 1.5 times larger 
than the sSFR for galaxies with M* > 1O 8 M over the 
redshift range z ~ 6 — 9 as shown in Fig. [3] (see also 
Sec. 3.4, iDaval fc Ferraral |2012| ) (note that galaxies with 
M* > 1O 8 M0 have not had time to evolve before z ~ 9 in 
the simulated volume). Examples of the evolution of sSFR 
for individual simulated galaxies of varying stellar mass are 
shown in Fig. IA3I in the appendix. 

More dramatic is the predicted trend with redshift 
above z = 6. As shown in Fig. [3l averaged over all LBGs 
with Muv < —15, the typical sSFR inferred from our sim- 
ulations rises by a factor ~ 4 from z ~ 6 to z ~ 12 (from 
sSFR =~ 4.5 Gyr" 1 to sSFR ~ 18.6 Gyr" 1 ). This is a 
strong prediction for future observations, and it is hearten- 
ing to note that, in the redshift range probed by current 
data (z ~ 6 — 7), our predicted sSFR values are in excel- 
lent agreement with the most recent observational results, 
the range of which is dominated by the degree of nebular 
emission-line correction applied. 



4 STELLAR MASS ASSEMBLY 

Having shown that the simulation successfully reproduces 
the key current observations of the galaxy population at 
z = 6 — 8, we can legitimately proceed to use our model to 
explore how the galaxies on the UV LF shown in Fig.[T]build 
up their stellar mass. We again use the galaxies at z ~ 6 as 
our reference point, and in what follows all galaxies with 
Muv < — 15 at z ~ 6 (which number 871) are colloquially 
referred to as z ~ 6 LBGs and their stellar mass is denoted 
as M*,z=6- 

Before discussing how z ~ 6 LBGs assemble their stel- 
lar mass, we briefly digress to explain how we identify their 
progenitors and the major branch of their merger trees be- 
tween z ~ 7 and z ~ 12, in snapshots spaced by Az ~ 1. We 
start our analysis from the simulation snapshot at z ~ 7. 
Since each star particle is associated with a bound struc- 
ture, in step one, by matching the star particles in the snap- 
shot at z ~ 7 and z ~ 6, we identify all the progenitors of 
each z ~ 6 LBG at z ~ 7; alternatively, each z ~ 7 galaxy 
has a z ~ 6 "descendant" galaxy associated with it. For 
each z ~ 6 LBG, the z ~ 7 progenitor with the largest to- 
tal (DM+gas+stellar) mass is then identified as the major 
branch of its merger tree at z ~ 7. In step two, matching the 
star particles in the simulation snapshots at z ~ 7 and z ~ 8, 
each z ~ 8 galaxy is associated to a "descendant" galaxy at 
z ~ 7 (i.e. we find all the progenitors of each z ~ 7 galaxy). 

Again, for each z ~ 7 galaxy, the z ~ 8 progenitor with 
the largest total mass is identified as the major branch of its 
merger tree at z ~ 8. Since we know the IDs of the major 
branch and all the progenitors of each z ~ 6 LBG at z ~ 7 
from step one above, and the ID of the z ~ 8 major branch 
and all progenitors of each z ~ 7 galaxy from step two, for 
each z ~ 6 LBG we can identify (a) all the progenitors at 
z ~ 7, 8, and (b) the major branch of the merger tree at 
z ~ 7, 8. This same procedure is then followed to find all the 
progenitors, and the major branch of the merger tree for 
each z ~ 6 LBG, at z ~ 9, 10, 11, 12; we end the calculations 
at z ~ 12 because only a few thousand star particles (and 



© 0000 RAS, MNRAS 000, 000-000 



Assembly of high-z galaxies 7 



"i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i/i r 


i i i i | i i i i | i i i i | i i 

- (b) z~8 / 

7 / K 

-- 

- /JSfe 

- / JHK V 

"i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 u r 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i 

- (c) z~9 

"i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i/i r 


- (d) z~10 / 

-- / x: 

/ *fc • • • 

i i i i i i i i i i i i i i i i i i i 


(e) z~ll - 
- / # - 

: / x: 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


- (f) Z~12 y/ 

: / . .r.r 

L . > 7 

i i i i l i i i i l i i i i l i i i i 



9 - 



8 - 



6 - 



7 8 9 

Log (M. z=6 /M ) 



Figure 4. The stellar mass of the major branch progenitor as a function of the z ~ 6 stellar mass, for z ~ 6 LBGs (small points) at the 
redshifts marked in each panel: z ~ 7 — 12. In each panel, the large symbols show the results for the 5 different example z ~ 6 LBGs 
discussed in the text: galaxy A (empty circle), B (empty square), C (empty star), D (asterix) and E (cross). The solid line indicates 
M*,MB = A4* iZ = 6 to guide the eye. 



a few hundred bound structures) exist at this epoch in the 
simulated volume, leading to very poor number statistics. 



4.1 Assembling the major branch 

We now return to our discussion on the stellar mass as- 
sembly of the major branch, M*,mb, for z ~ 6 LBGs. As 
expected from the hierarchical structure formation scenario 
where larger systems build up from the merger of smaller 
ones, the earlier a system starts assembling, the larger its 
final mass is likely to be; alternatively, this implies that the 
progenitors of the largest systems start assembling first, with 
the progenitors of smaller systems assembling later. A natu- 
ral consequence of this behaviour is that it leads to a corre- 
lation between M» jZ= 6 and M* t MB, i.e. at any given z, the 
major branch stellar mass is generally expected to scale with 
the final stellar mass at z ~ 6. 

This behaviour can be seen clearly in Fig. [4] firstly, we 
see that the progenitors of the most massive z ~ 6 LBGs 
start assembling first, and progenitors of increasingly smaller 
systems appear with decreasing z. Indeed, while the pro- 
genitors of z ~ 6 LBGs with M* ^ 1O 8 M0 already exist at 
2 ~ 12, there is a dearth of progenitors of the lowest-mass 



galaxies, with M* ~ 1O 7 M0, which start building up as late 
as z ~ 9. It can also be seen from the same figure that there 
is a broad trend for the z ~ 6 LBGs with the largest M*, z =6 
values to have largest values of M t| MS at all the redshifts 
studied z ~ 7— 12, although the scatter grows substantially, 
and only a few hundred progenitors exist at z > 11. The 
wide spread in M*,mb for a given final M» jZ= 6 value, points 
to the varied stellar mass build-up histories of the galaxies 
in the simulation; the relation between M*,mb and M* !Z =6 
inevitably tightens with decreasing redshift (see also Table 

E>. 

Finally, we note that, averaged over all z ~ 6 LBGs, 
only about 0.3% of the final stellar mass has been built up 
by z ~ 12. The major branch of the merger tree steadily 
builds up in mass at a rate (mass gained/Myr) that increases 
with decreasing redshift such that about (1, 3, 8, 19, 37)% of 
the final stellar mass is built up by z ~ (11,10,9,8,7) as 
shown in Table 1; this implies that z ~ 6 LBGs gain the 
bulk of their stellar mass (« 63%) in the ~ 150 Myr be- 
tween z ~ 7 and z ~ 6, and only about a third is assembled 
in the preceding ~ 400 Myr between z ~ 12 and z ~ 7. This 
rapidly increasing stellar mass growth rate with decreasing 
z probably results from negative mechanical feedback be- 
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Figure 5. The ratio of the stellar mass of the major branch progenitor and all the other progenitors (i.e. excluding the major branch 
progenitor) for each z ~ 6 LBG, as a function of the z ~ 6 stellar mass. The panels show the results for different redshifts, z ~ 7— 12, and 
the small grey (green) points show the results for LBGs with more than one (a single) progenitor. Since galaxies with a single progenitor 
have an undefined value of M. t MB /^*,allprogy we use an arbitrary ratio of 10 2 to show these galaxies. The large symbols show the 
results for the 5 different example z ~ 6 LBGs discussed in the text: A (empty circle), B (empty square), C (empty star), D (asterix) 
and E (cross). 



coming less important as galaxies grow more massive: even 
a small amount of star formation in the tiny potential wells 
of early progenitors can lead to a partial blowout/full blow- 
away of the interstellar medium (ISM) gas, suppressing fur- 
ther star formation (at least temporarily). These progenitors 
must then wait for enough gas to be accreted, either from 
the IGM or through mergers, to re-ignite star formation; 
negative feedback is weaker in more massive systems due to 
their larger potential wells (see also lMaio et al.ll201lT ). 

To help illustrate and clarify the above points, we show 
the stellar mass growth of five different z ~ 6 LBGs with 
log[M*,z =6 /M ] = (6.7,7.6,8.3,9.1,9.8); these are labelled 
(A,B,C,D,E) respectively. As seen from panel (f) of Fig. [4] 
the progenitors of the two most massive example galaxies, 
D and E, appear as early as z ~ 12, making some of their 
stellar content at least as old as 550 Myr by z ~ 6. How- 
ever, although these galaxies started building up their stellar 
mass early, their stellar mass-weighted ages are, of course, 
substantially younger, at (227, 132) Myr respectively; con- 
sistent with the statement made above, the bulk of their 
stellar mass was formed in the 150-200 Myr immediately 



prior to z ~ 6 and these galaxies have only assembled about 
(1.9, 0.1)% of their final stellar mass at z ~ 12. The value of 
M- Mi mb steadily builds up, either by merging with the other 
progenitors of the "descendant" galaxy, or due to star for- 
mation in the major branch itself such that, by z ~ 10, these 
galaxies have increased in mass by a factor of about 4 (6). At 
z ~ 10 (9), the progenitor of galaxy C (B) appears, putting 
the age of the oldest stars in this galaxy at > 450 (380) Myr 
by z ~ 6; again the bulk of the stellar mass for both these 
galaxies forms in the ~ 156 Myr between z ~ 7 and z ~ 6. 
The progenitor of galaxy A finally appears at z ~ 7, having 
formed a tiny stellar mass of M*,mb — 1O 6 M0. Although 
by z ~ 7, the five galaxies illustrated have individually built 
up between 15-50% of their final stellar mass, their average 
z ~ 7 stellar mass is about 30% of the final value; this is 
very consistent with the global average value of 37%. 
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Figure 6. The stellar mass functions for the progenitors of z ~ 6 LBGs: the left and right-hand panels show the results for the major 
branch of the merger tree, and all the other progenitors (i.e. excluding the major branch), respectively. In each panel, the lines show the 
results at different redshifts: z ~ 7 (dotted red), z ~ 8 (short dashed blue), z ~ 9 (long dashed green), z ~ 10 (dot-short dashed orange), 
z ~ 11 (dot-long dashed cyan), 2 ~ 12 (short dashed- long dashed violet). For comparison, the solid line in each panel shows the stellar 
mass function for z ~ 6 LB Gs. In the left-hand pa nel, filled (empty) points show the corrected (uncorrected) z ~ 6 stellar mass functions 
inferred observationally bv lGonzalez et all i201ll ). 



4.2 Relative growth: major branch versus all 
progenitors 

We now ask how the growth of M»,mb compares to the 
stellar mass growth of all the other progenitors (i.e. exclud- 
ing the major branch itself) of z ~ 6 LBGs, M t ^ a u prog . 
Naively it might be expected that when the major branch 
of a galaxy first forms, its stellar mass might be low com- 
pared to that contained in all the other progenitors; the 
ratio Mt^MB/Mt^aiiprog would increase with time as M»,mb 
builds up, becoming the dominant stellar mass repository. 
This is indeed the general behaviour found for the progen- 
itors of z ~ 6 LBGs, where the value of the total stellar 
mass contained in all the major branch progenitors of z ~ 6 
LBGs, EM*, ms, increases from 1.8 to 6.1 times the total 
mass contained in all the other progenitors, T>M, : aii P rog, as 
the redshift decreases from z ~ 12 to z ~ 7 (as shown in 
Fig. [S] see also Table [TJ; while between z ~ 12 and z ~ 10, 
all the other progenitors put together hold almost half as 
much stellar mass as the major branch, this value rapidly 
decreases thereafter, such that by z ~ 7, the other progen- 
itors together contain only about 15% of the mass in the 
major branch. 

The ratio M 4] mb /M r ^ a u prog increases very slowly be- 
tween z ~ 12 and z ~ 10 because, at these redshifts, both 
the major branch and all the other progenitors seem to be 
growing in stellar mass at the same pace. However, at z ^ 10, 
the major branch progenitor starts growing more quickly, 
possibly due to the fact that, as the major branch becomes 
increasingly more massive, it is increasingly less affected by 
the negative mechanical feedback (i.e. SN powered outflows) 



that can quench star formation in the other smaller progen- 
itors. 

Further, as discussed above in Sec. 14.11 the progeni- 
tors of successively less massive z ~ 6 LBGs appear with 
decreasing z, leading to a scatter of more than an or- 
der of magnitude in M*,a/s / 'M, >aUprog , at any of the red- 
shifts shown in Fig. [S] due to their varied assembly histo- 
ries. Indeed, although on average the major branch is the 
heavyweight in terms of the stellar mass content, at any 
of the redshifts studied, there are always a few galaxies 
where the major branch only contains about 10% of the 
total stellar mass at that redshift. Finally, we note that, 
at all redshifts, a number of progenitors show a value of 
M* t MB = M* t a.ii P rog'- these are galaxies with a single pro- 
genitor, the major-branch one. Although such galaxies do 
not have a defined M,,mb /M* ia llprog value, we have used 
an arbitrary value of 10 2 to represent them in Fig. \5\ 

The five example LBGs discussed in Sec. 14.11 help to 
illustrate these points: galaxy E, the most massive galaxy 
in the simulation box, assembles from tens of tiny progen- 
itors, and at z ~ 12 its major branch contains only about 
25% of the stellar mass in place at that epoch, with the 
bulk (75%) contained in all the other progenitors. On the 
other hand, galaxy D has a very different assembly history; 
at z ~ 12 it already has major branch progenitor that is 
about three times more massive than all the other progen- 
itors put together. Between z ~ 12 and z ~ 9 the value 
of M tt MB/M* t aiiprog remains almost unchanged for E; al- 
though M*,mb has doubled for E between these redshifts. 
This is an example of a case where all the other progen- 
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itors put together grow in mass as rapidly as the major 
branch. For galaxy C, on the other hand, all the other pro- 
genitors grow more quickly than the major branch between 
z ~ 10 and z ~ 8, leading to a decrease in the value of 
M,,mb /M* :a iip r0 g with time. Galaxy B is a classic exam- 
ple of a galaxy that starts off with a single progenitor (the 
major branch one) at z ~ 9, and then briefly has two progen- 
itors at z ~ 8, before these merge so that the major branch 
again holds all the stellar mass by z ~ 7. At z ~ 7, the re- 
cently formed major branch of A already holds the majority 
of the total stellar mass assembled; by this redshift, D and 
E have assembled a major branch that is about 30 and 6 
times heavier in stellar mass than all the other progenitors 
together. 

We now summarize the relationship between M* tZ =B t 
M* t MB and M+ >a uprog by representing these quantities in 
terms of their mass functions at all the redshifts studied (see 
Fig. H3). We start by noting that, while the theoretical z ~ 6 
LBG stellar mass function is in broad agreement with the 
observed one for M* > 1O 9 M0, it is much steeper than cur- 
rently inferred from the observations for lower values of M t . 
This may be due to the fact that the observational stellar 
mass function has been constructed by inferri ng M* from the 
UV l uminosities of LBGS with Muv < —18 (Gonzalez et al. 
Son]), and is incomplete at the low-mass end. Further, the 
theoretical stellar mass function of z ~ 6 LBGs peaks at 
M* ~ 10 7 ~ r,5 M Q , as shown in panel (a) of Fig. g] The de- 
crease in the number density of the more massive objects 
is only to be expected from the hierarchical structure for- 
mation model. The decrease in the number density of lower 
stellar mass systems, on the other hand, arises from the fact 
that due to their low stellar masses, not many of these galax- 
ies produce enough UV luminosity to be visible as LBGs 
with Muv < —15, i.e. this drop is not an artefact of the 
simulation resolution; the requirement that a galaxy has at 
least 10 star particles corresponds to a resolved stellar mass 
of 10 6 M Q . 

In panel (a) of Fig. HJ] we show the evolving stellar mass 
function of the major branch progenitors at z ~ 7 — 12 
which, as expected, shifts to progressively lower M* values 
with increasing redshift. The major branch progenitors com- 
pletely dominate the high mass end of the evolving stel- 
lar mass function, as becomes clear by comparison with 
panel (b) which shows the stellar mass functions of all the 
other progenitors (i.e. excluding the major branch) of z ~ 6 
LBGs. This shows that (i) these are always less massive than 
10 8,5 Mq at any of the redshifts considered, with this mass 
threshold decreasing with increasing z, as expected, and (ii) 
the contribution of the other progenitors is largest at the 
highest z ~ 10 — 12 where these contain about half the stel- 
lar mass in the major branch. When compared, it is clear 
that the major branch dominates over the contribution from 
all the other progenitors for M» ^ 10 7 Mq. However, at the 
low mass end (M. < 10 7 M Q ), there is an abundance of tiny 
progenitors that contributes to the mass function by a factor 
of about 1.25 (0.7) compared to the major branch for z ~ 12 



4.3 Mass growth driver: star formation or 
mergers? 

We now return to Figs. [4] and [5] to consider an interest- 
ing question: what is the relative contribution of merg- 
ers/accretion (Mace) versus star formation (M s f) to build- 
ing up the total stellar mass of the major branch? A first 
hint of the answer can be obtained by comparing the pop- 
ulation average numbers as presented in columns 1 and 3 
of Table 1 (see also Sees. O and g2J) for z ~ 6 LBGs: at 
z ~ 7, the total stellar mass in the major branch has a value 
0.37M» jZ= 6, and all the other progenitors contain a total 
stellar mass of 0.06M» iZ =6- This implies that the merger of 
all the other progenitors into the major branch can increase 
the mass up to 0.43M* ;Z= 6, and thus that the remaining 
57% of M ti 2 = 6 must be produced by star formation in the 
redshift interval z ~ 7 to z ~ 6. This then implies that, 
on average, in the redshift interval z ~ 7 — 6, stellar mass 
growth due to star formation dominates growth due to the 
accretion of pre-existing stars in the other progenitors by a 
factor 0.57/0.06 ~ 10. 

For a more complete answer, we have carried out this 
calculation on a galaxy- by-galaxy basis. At any redshift, the 
stellar mass accreted by the major branch due to mergers 
between redshifts z and z — 1 is calculated as the differ- 
ence between the mass of all the major branch progeni- 
tors and the major branch mass at z, i.e. M aC c(z,z — 1) = 
M, a i iprogMB(z) — Af* ; mb (z) ■ The stellar mass assembled 
due to star formation is calculated as the difference between 
the major branch mass at 2 — 1, and the mass of all the ma- 
jor branch progenitors (including the major branch itself) 
at z, i.e. M af (z,z — 1) = M*,mb(z - 1) - M* ia llprogMB(z), 
where M» ta ii prog MB is the stellar mass in all the progenitors 
of the major branch. 

We note here that while the 'mass growth due to star 
formation' M s f so calculated is an accurate estimate of the 
amount of new stellar mass contributed to the final z ~ 6 
galaxy in the preceding Az = 1, we cannot be sure that all 
this star-formation takes place 'within' the major branch. 
Some of this star formation could clearly take place in the 
minor progenitors before they finally merge with the ma- 
jor branch at some point between z and z — 1. However, 
to better define what fraction of this star-formation activity 
takes place within the major branch would require the anal- 
ysis of simulation snapshots very closely spaced in z; since 
such snapshots have not been stored, this calculation would 
require re-running the simulation which is not feasible at 
present. In any case, regardless of the precise location of the 
star-formation activity, this calculation delivers a robust es- 
timate of the contribution of star formation in the preceding 
A 2 = 1 to the stellar mass of the galaxy. 

The main results of these calculations are as follows. 
In all redshift intervals, for the majority of galaxies, we 
find most of the mass growth results from star formation 
in the major branch, with mergers bringing in (relatively) 
tiny amounts of stellar mass such that M s f » M acc as can 
be seen from Fig. [7]and Table 1; indeed, as is seen from the 
same table, ~ 90% of the total stellar mass of z ~ 6 LBGs 
is assembled by star formation in the major branch, with 
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Figure 7. The fraction of the total stellar mass assembled in z ~ 6 LBGs by star formation in the major branch (grey circles), and by 
merging of the minor progenitors (green squares), in the redshift ranges marked in each panel, as a function of the final stellar mass of 
each object at z ~ 6. The large squares (circles) in each panel show the mass assembled by the 5 example z ~ 6 LBGs by mergers (star 
formation) in each redshift interval: A (purple), B (black), C (maroon), D (blue) and E (red). Note that when any galaxy/major branch 
progenitor has only one progenitor at the previous redshift, all the mass is assembled because of star formation; the mass brought in by 
mergers is, by definition, zero. 



mergers contributing only ~ 10% to this final value. Never- 
theless, while this average behaviour is clear, there is a wide 
variety in the history of individual galaxies, and at almost 
all 2, several galaxies in fact show values of M acc « M s f. 
Finally, our result that the bulk of M*, z= 6 has to be built in 
the ~ 150 Myr between z ~ 7 and z ~ 6 is verified by panel 
(a) of Fig.[7]where star formation in the major branch leads 
to an average of about 56% of M* iZ= 6 being formed between 
these redshifts. 

We can try to further clarify these results by again dis- 
cussing the histories of our five example galaxies: as shown 
in Fig. [4] galaxies D, E have built up only about 3, 0.3% of 
M* iZ =6 by z ~ 11; from panel (f) of Fig. [7] we can clearly 
see that while for E this mass has solely been built up by 
M s f, for D, there is some significant contribution due to 
Macc- These galaxies then slowly build up their stellar mass, 
mostly through star formation until z~9 — 8, 8 — 7 when 
mergers of smaller systems start making a noticeable con- 
tribution; mergers contribute as much to the stellar mass 
as star formation in the major branch for C, D and E at 
z ~ 8 — 6, 8 — 7 and 9 — 8, respectively. On the other hand, 



galaxy E only has a single progenitor for most of its life and 
thus grows solely by star formation in the major branch, ex- 
cept between z ~ 8 and z ~ 7 when it has two progenitors 
(see also fig. [5} and merger of the other major-branch pro- 
genitor contributes a small amount to M*, z =6- Finally, about 
40-70% of M*, z= 6 for these five example galaxies is built-up 
by star formation in the ~ 150 Myr between z ~ 7 and z ~ 6, 
with mergers contributing only about 2-25% to A/* jZ= 6 (the 
remaining mass pre-existing in the major branch at z — 7). 



5 UV LUMINOSITY EVOLUTION 

Now that we have discussed how z ~ 6 LBGs assemble their 
stellar mass, we can also study how their UV luminosity 
varies over time. We address this question by comparing 
the absolute UV magnitude of the major branch progeni- 
tor (Muv,mb) at each redshift, with the final absolute UV 
magnitude of the z ~ 6 LBGs, Mjyv>=6, as shown in Fig. [8] 
Although (E(B — V)} < 0.05 at z ~ 6 and is expected to de- 
crease further with increasing z as a result of galaxies being 
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Figure 8. The absolute UV magnitude of the major branch progenitors with Muv, MB < — 15 for each z ~ 6 LBG (small points), as a 
function of the final absolute UV magnitude at z ~ 6. The panels show the results for different redshifts, z ~ 7 — 12, as marked, and the 
large symbols show the results for the example z ~ 6 LBGs discussed in the text: B (empty square), C (empty star), D (asterix) and E 
(cross); galaxy A does not appear on this plot since its progenitors are always fainter than the limiting value of Muv < — 15 shown here. 



younger and less metal enriched, we have nonetheless calcu- 
lated the dust enrichment of each galaxy in all the snapshots 
used between z ~ 12 and z ~ 6, consistently applying the 
dust model described in Sec. 12.11 

The first result we find is that Muv.mb oc Muv,z=6 
at least for z ~ 7 — 10, consistent with the most massive 
z ~ 6 LBGs having the most massive major branch pro- 
genitor; this correlation becomes sketchy at z 11 as a 
result of the scarcity of progenitors brighter than the ap- 
plied magnitude limit of Muv < — 15 above this redshift. 
The average absolute magnitude shifts towards lower val- 
ues (i.e. increasing luminosity) with decreasing z, thus pro- 
ducing the effect of luminosity evolution in the luminosity 
function (see below). Moreover, the brightest sources appear 
to display fairly steady luminosity evolution on a source- 
by-source basis; for example, the progenitors of the z ~ 6 
galaxies with Muv < —18.5 constitute ~ 90% of the galax- 
ies with Muv < — 18 at z ~ 7, and still provide ~ 75% of 
the galaxies brighter than Muv — —18 at z ~ 8. 

However, at fainter luminosities the situation is more 
complex, as illustrated by the fact that between z ~ 7 and 
z ~ 10 a number of individual major-branch progenitors 
with M* jZ= 6 < 1O 8 M0 are more luminous than their (higher 



mass) z ~ 6 descendants (i.e. Muv,mb < Muv,z=6, as shown 
in panels (a) to (d) of Fig.|HJ. This result may seem puzzling, 
given that even at z ~ 7, these small galaxies have only been 
able to assemble a stellar mass M*,mb ~ 0AM* )Z= e, but is 
explained by the more stochastic star-formation histories of 
the lower-mass galaxies (as illustrated in Fig. IA2l in the ap- 
pendix). In other words, a short-lived burst of star- formation 
in a low mass major-branch progenitor, which decays on 
times-scales of a few Myr, can easily produce a temporary 
enhancement in UV luminosity which exceeds that of the 
(still fairly low mass) descendant at z ~ 6. 

As for the 5 example galaxies that have been discussed 
above, the smallest at z ~ 6, galaxy A is right at the limit of 
our magnitude cut, with Muv,z=6 ~ —15.1, and its major- 
branch progenitors are never luminous enough to clear this 
threshold at higher redshift. The major branches of B and 
C (Muv,z=6 — —16.6, —17.9) already have a stellar mass 
M*,mb — (15, 6)%M« jS= 6 respectively, when they become 
visible with Muv.mb — —15.2 at z ~ 8. For D and E, which 
have final z ~ 6 magnitudes of Muv,z=e — —18.8, —21 re- 
spectively, the major-branch progenitors are visible at all the 
redshifts z ~ 7 — 12 with the UV magnitude decreasing (i.e. 
their luminosities brightening) monotonically with decreas- 



© 0000 RAS, MNRAS 000, 000-000 



Assembly of high-z galaxies 13 




Figure 9. UV LFs for the major branch of the merger tree for z ~ 6 LBGs, for z ~ 7 to z ~ 12, as marked in each panel. In each 
panel, the long dashed (red), short dashed (green) and dotted (blue) lines show the major branch of the merger tree for z ~ 6 LBGs 
with Muvz=6 < — 15, Muv z= 6 < —18 and —18 ^ Muv,z=6 ^ — 15, respectively. For comparison, the solid (black) line in each panel 
shows the UV LF obtained using all the galaxies in the lOcMpc simulation box with Myy < —15, at the redshift marked in that panel. 



ing z. However, we note that the UV luminosity of E grows 
faster compared to D for z ^ 9 as a result of its faster stel- 
lar mass growth that naturally produces more UV photons, 
consistent with our finding above that most mass growth in 
a given redshift interval is driven by star formation. 

Finally, we show, and deconstruct the simulated UV 
LFs for LBGs in the simulation with Muv < — 15, over the 
redshift range z ~ 7—12. The full simulated UV LFs at 
z ~ 6, 7, 8 have already been shown in Fig. [TJ but in Fig. 
[9] we extend this to predict the form of the total UV LF 
up to z ~ 12, and also separate out the contribution of 
the major-branch progenitors to the LF. At all the redshifts 
shown (z ~ 7 — 12), the major branch progenitors of the 
z ~ 6 LBGs dominate the UV LF over the entire magnitude 
range probed here. At the brightest magnitudes they are, 
not surprisingly, completely dominant, and even at the faint 
end (Muv ~ — 16) they constitute ~ 60% of the number 
density of objects in the LF. 

To further clarify the physical evolution of the galaxy 
population, we have also produced the UV LF for the major- 
branch progenitors of z ~ 6 LBGs with (a) Muv,z=6 < — 18 
and (b) —18 ^ Muv,z=e ^ —15. We find that the major- 



branch progenitors of the former dominate the high lumi- 
nosity end of the UV LF at all z and make a negligible 
contribution to the faint end. By contrast, the progenitors 
of the faintest z ~ 6 LBGs with — 18 ^ Muv,z=6 ^ — 15 
contribute to the faint end of the UV LF at all the redshifts 
studied. 



In conclusion, our simulation suggests that it is indeed 
reasonable to expect a steady brightening of the bright end 
of UV LF during the first billion years, and that this is pri- 
marily driven by genuine physical luminosity evolution (i.e. 
steady brightening, albeit not exponential) of a fixed subset 
of the highest-mass LBGs. However, at fainter magnitudes 
the situation is clearly much more complex, involving a mix 
of positive and negative luminosity evolution (as low-mass 
galaxies temporarily brighten then fade) coupled with both 
positive and negative density evolution (as new low-mass 
galaxies form, and other low-mass galaxies are consumed by 
merging) . 
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6 CONCLUSIONS 

We have used state-of-the-art high-resolution SPH cosmo- 
logical simulations, specially crafted to include the physics 
most relevant to galaxy formation (star formation, gas cool- 
ing and feedback) and a new treatment for metal enrichment 
and its dispersion that allows us to study the transition from 
metal free PopIII to metal enriched PopII star formation, 
and hence simulate the emergence of the first generations of 
galaxies in the high-redshift universe. By combining simu- 
lation snaps hots at z — 6 — 1 2 with a previously developed 
dust model (|Daval et al.l2lll<:i. in addition to calculating the 
LBG UV LFs, SMD and sSFR for the current magnitude 
limit of Muv — —18, we have extended our results down to 
Muv — —15 in order to make specific predictions for up- 
coming instruments such as the JWST. The main results 
from our comparison of the population predictions from the 
theoretical model and the constraints provided by current 
data can be summarized as follows. 

(i) UV LFs: we find that the simulated LBG UV LFs 
match both the amplitude and the slope of the observations 
in the current range of overlap; the faint-end slope is found 
to remain almost constant between z ~ 6 and z ~ 8, with a 
value consistent with a ~ —2. 

(ii) SMD: the SMD for the simulated LBGs with M uv < 
— 18 is in extremely good agreement with the observed values 
at z ~ 6 — 8. Interestingly, we find that as a result of their 
much larger number density, about 1.5 times the currently- 
observed stellar mass is locked up in the undetected, faint 
LBGs with -18 < Muv < -15 at these redshifts. 

(iii) sSFR: we find that the sSFR decreases gently with 
increasing M, values since even a small amount of star for- 
mation in galaxies with low M* values is enough to push 
up the sSFR. The theoretical sSFR increases from about 
4.5 Gyr -1 at z ~ 6 to about 11 Gyr -1 at z ~ 9. In the cur- 
rently accessible redshift range of overlap (z ~ 6 — 7), the 
simulated sSFRs are in excellent agreement with those in- 
ferred from the latest observations (especially after correct- 
ing for emission- line contributions), but the above figures 
obviously provide a strong prediction of rapidly increasing 
typical values of sSFR at higher redshifts, a prediction po- 
tentially testable with the JWST. 

Having demonstrated that the population predictions 
from the simulations match well with the existing observa- 
tional data, we have validated our model as a potentially 
viable description of the physical growth and evolution of 
high-redshift galaxies. We have therefore proceeded to ex- 
amine the detailed galaxy-by-galaxy behaviour within the 
simulation in an attempt to gain physical insight into the 
observed population statistics. To do this we have identi- 
fied all the progenitors as well as the major branch of the 
merger tree for z ~ 6 LBGs (Muv < —15) at higher red- 
shift z ~ 7, 8, 9, 10, 11, 12. This has enabled us to study how 
these z ~ 6 LBGs have assembled their stellar mass over the 
preceding ~ 600 Myr of cosmic history, and also investigate 
how their UV luminosity has varied during this time. 

The story that emerges regarding the lives of these 
galaxies can now be summarized as follows. At z ~ 12, the 



progenitors of z ~ 6 LBGs contain, on average, only 0.3% of 
the final z ~ 6 stellar mass and the rate of stellar- mass build- 
up increases with decreasing redshift such that by z ~ 7, 
these progenitors contain about 37% of the total stellar mass 
contained in z — 6 LBGs; this increasing SF efficiency ap- 
pears to be the result of the decreasing effects of negative 
feedback (blowout/blowaway) as the potential well of the 
major branch grows. These numbers imply that the bulk 
(~ 56%) of M, is formed in the ~ 150 Myr between z ~ 7 
and z ~ 6. Further, we find, unsurprisingly, that the pro- 
genitors of the most massive z ~ 6 LBGs start assembling 
first; while progenitors of z ~ 6 LBGs with M* > 1O 8 M0 
exist as early as z ~ 12, progenitors of the lowest mass 
(M» 1O 7 M0) z ~ 6 galaxies appear as late as z ~ 9. As a 
result of their early assembly, the most massive z ~ 6 LBGs 
generally have the most massive major branch progenitor at 
the redshifts studied. 

Compared to all the other progenitors (i.e. excluding 
the major branch itself), the major branch is always the 
dominant M* component; at z — 12 the major branch typ- 
ically contains about twice the total mass across all other 
progenitors of z ~ 6 LBGs, and by z ~ 7 this ratio has risen 
to about six. In the redshift range z ~ 12 — 10, the major 
branch and all the other progenitors grow at approximately 
the same rate, but thereafter the major branch grows more 
quickly in M* , possibly due to the decreasing impact of neg- 
ative feedback. However, while the global trend is clear, at 
every redshift there are outliers in which the major branch 
contains as little as 10% of the total stellar mass assem- 
bled by that redshift, highlighting the varied histories of the 
galaxies in the simulation. 

We have also attempted to determine the relative im- 
portance of mass added by pre-existing stars in merging pro- 
gentors, as compared to that added by new star-formation 
activity in a given redshift interval. To compute this we have 
defined the stellar mass 'accreted' by the major branch at 
redshift z — 1 as the stellar mass of all its other progenitors 
at 2 and the mass gained by 'star formation' as the differ- 
ence of the stellar mass of the major branch at z — 1 and all 
its progenitors at z. From this calculation, we have shown 
that ~ 90% of the total z ~ 6 stellar mass is built up by star 
formation, and mergers contribute only ~ 10% to the final 
stellar mass. However, the average contribution of mergers 
increases with decreasing redshift, from ~ 0.01% at z ~ 12 
to ~ 6% at z ~ 6, as more and more progenitors form and 
fall into the ever-increasing major-branch potential well. 

Finally, we have tracked the variation of UV luminosity 
in the galaxies and their various progenitors over the red- 
shift range z ~ 12 — 6. At least since z ~ 9, we find that it 
is largely the same set of relatively massive galaxies which 
occupy the bright end of the UV luminosity function, and 
that these objects brighten steadily with time (albeit not ex- 
ponentially), thus providing a physical basis for interpreting 
the evolution of the bright end of the LF in terms of lumi- 
nosity evolution. However, at fainter luminosities the situa- 
tion is more complex, as illustrated by the fact that between 
z ~ 7 and z ~ 10 a number of individual major-branch pro- 
genitors with M*, z= 6 < 1O 8 M0 are more luminous than their 
(higher mass) z ~ 6 descendants; this appears to be one con- 
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sequence of the more stochastic star-formation histories of 
the lower-mass galaxies in our simulation (as illustrated by 
the examples of individual galaxy formation histories pro- 
vided in the appendix). 

In conclusion, our simulation suggests that it is indeed 
reasonable to expect a steady brightening of the bright end 
of UV LF during the first billion years, and that this is 
primarily driven by genuine physical luminosity evolution 
(i.e. steady brightening) of a fixed subset of the highest- 
mass LBGs. However, at fainter magnitudes the situation 
is clearly much more complex, involving a mix of posi- 
tive and negative luminosity evolution (as low-mass galaxies 
temporarily brighten then fade) coupled with both positive 
and negative density evolution (as new low-mass galaxies 
form, and other low-mass galaxies are consumed by merg- 
ing). Thus, it is probably not reasonable to aim to explain 
the entire evolution of the UV LF, over a broad range in 
absolute magnitude, in terms of either pure luminosity evo- 
lution or pure density evolution. However, improved data 
are clearly required to distinguish between the viability of 
the model presented here, and alternative scenarios. 

We end by discussing a few caveats: first, the simula- 
tion has been run using a star formation density threshold 
of 70 cm -3 . Obviously, the higher this threshold, the later is 
the onset of star formation as the gas needs more time to 
condense. However, exploring density t hresholds ranging be- 
tween 0.1 — 70 cm~ 3 . lMaio et ah! (|2010t ) have shown that the 
the SFRs obtained from all these runs converge as early as 
z ~ 13. Second, ejection of both gas and metal particles into 
the IGM remains a poorly understood astrophysical process. 
Although in our simulations, winds originate as a result of 
the kinetic feedback from stars, scenarios like gas stripping, 
shocks and thermal heating from stellar processes could play 
an important role, given that it is extremely easy to lose 
baryonic content from the small DM potential wells of high- 
z galaxies. It is hoped that in the future, an improvement 
in our theoretical understanding, and further inputs from 
observations will shed further light on these issues. 
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APPENDIX A: EXAMPLES OF STELLAR 
MASS ASSEMBLY, SFR AND SSFR FOR Z ~ 6 
LBGS 

We show the z-dependent stellar mass growth, and the 
SFR and sSFR evolution for the major branch of the 
merger tree for 18 'typical' LBGs selected from the simu- 
lation at z ~ 6, with M*, z =e ranging between 10 7 '°" 9 - 8 M Q . 
In these illustrative plots, the SFR has been calculated as 
M t = [M,(t) - M,(t - 1)]/At where M*(t) - M*(t - 1) 
shows the stellar mass built within a certain time step and 
At is the width of the time-step. Note that this stellar mass 
includes both the mass built by star formation inside the 
major branch, as well as the mass gained by merging with 
progenitors of the major branch (see sec. 14. 3[) . 
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Figure Al. The total stellar mass assembled in the major branch by the redshift shown on the x-axis, either through mergers or star 
formation. Each of the panels shows the cumulative mass for different galaxies, with the M t}Z —% value marked. 
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Figure A2. The SFR computed as the ratio of the stellar mass assembled by the major branch in a given time step, by the time-step 
width, as a function of the redshift. Each of the panels shows the SFR for different galaxies with the final M* >z — g value marked. 
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Figure A3. The sSFR computed as the ratio of the SFR to the stellar mass assembled by the major branch, as a function of the 
redshift. Each of the panels shows the sSFR for different galaxies w ith the final M*, z =6 v alue marked. The horizontal thick (thin) lines 
show the average (spread) of the sSFR computed observationally dGonzalez et al.ll2010l) between z ~ 4 and z ~ 8; to guide the eye, 
these observational results have been extended to z = 16 assuming that the observed sSFR remains constant to such redshifts. 
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